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Abstract 



We present precision corrections to dispersion relation bounds on form 



^ . factors in bottom hadron semileptonic decays and analyze their effects on 

H ! 

parameterizations derived from these bounds. We incorporate QCD twodoop 
and nonperturbative corrections to the two-point correlator, consider form 
factors whose contribution to decay rates is suppressed by lepton mass, and 
implement more realistic estimates of truncation errors associated with the 
parameterizations. We include higher resonances in the hadronic sum that, 
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together with heavy quark symmetry relations near zero recoil, further tighten 
the sum rule bounds. Utilizing all these improvements, we show that each 
of the six form factors in B — ► Dlv and B — ► D*lv can be described with 
3% or smaller precision using only the overall normalization and one unknown 
parameter. A similar one-coefficient parameterization of one of the A& — > A c lv 
form factors, together with heavy quark symmetry relations valid to order 
1/m 2 , describes the differential baryon decay rate in terms of one unknown 
parameter and the phenomenologically interesting quantity Aa Ma 6 — mj,. 
We discuss the validity of slope-curvature relations derived by Caprini and 
Neubert, and present weaker, corrected relations. Finally, we present sample 
fits of current experimental B — »■ D*lv and B — >■ Dlv data to the improved 
one-parameter expansion. 
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Typeset using REVTpjX 



I. INTRODUCTION 



The decays of bottom hadrons provide fertile ground for explorations of both weak and 
strong interactions. The empirical smallness of the CKM elements \V c b\ and \V u b\ implies 
that 6-hadron decays are relatively long-lived, and thus more easily analyzed, while the 
heaviness of the b quark means that the heavy quark effective theory (HQET) JT|J^] provides 
a reliable expansion for describing the approximate decoupling of the b quark from the rest 
of the hadron, leading to a substantial simplification of the strong-interaction dynamics. 

Semileptonic decays of 6-hadrons are particularly tractable from the theoretical point 
of view, since the leptonic current may be trivially separated from the hadronic current. 
Furthermore, focusing upon exclusive decays permits one to avoid questions regarding the 
validity of quark-hadron duality near kinematic endpoint regions. On the other hand, ex- 
clusive processes must be described in terms of a number of nonperturbative form factors 
that encapsulate the physics of the hadronization process. 

HQET has provided a substantial leap forward in the exclusive analysis, demonstrating 
that heavy-to-heavy quark transitions at zero recoil are accompanied by a complete overlap of 
initial- and final-state hadron wave functions • Consequently, each form factor possesses 
a well-defined normalization at this kinematic point. For the case of — > £)W transitions, 
for example, each of the 20 form factors either vanishes or is proportional to a single universal 
form factor £, which represents this wave function overlap as a function of momentum 
transfer. Sum rule constraints restrict the slope and magnitude of these form factors £§ . 
but otherwise shapes of the form factors are a priori undetermined functions of momentum 
transfer. 

However, some recent work [0-18] has demonstrated that one may obtain rather nontriv- 
ial and model-independent constraints on the shape of such form factors in the context of 
dispersion relations. Using basic field theoretic properties like causality and crossing sym- 
metry, one finds that the shape of a given form factor is determined to high accuracy over 
its entire kinematic range by its value at only a small number of points [fZ| P~Q|| - Indeed, for 
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the case of B — > D^Hu, the HQET normalization plus only one additional point determines 
most of the form factors to 3-7% accuracy |7]]. Similar analyses provide interesting and 



useful constraints on B -> wiu 0,[lO|-[T3l, B -> p£z/ 0], and A fe -> A c £z/ 0. 

At the core of this analysis lies a dispersion relation in momentum space that relates the 
integral of an inclusive production rate to a two-point function evaluated far from physical 
resonances, where perturbative QCD is quantitatively reliable. The production rate is a 
sum of positive-definite exclusive rates; it follows that all contributions to the hadronic side 
serve to saturate more and more closely the bound from the perturbative side. Clearly, two 
improvements to this procedure include a better calculation of the QCD side (which includes 
both multi-loop effects and nonperturbative vacuum condensates), and the inclusion of more 
states on the hadronic side. The former may be extracted from existing literature, while 
the latter may be achieved using heavy quark symmetry relations. The inclusion of higher 
states has previously been investigated for elastic B — ► B scattering [T/J and for B — > -niv 
decays [|ll|], but not for b — > c transitions. For mesonic and baryonic b —>■ c transitions, 
the inclusion of higher states improves considerably the precision of the dispersion-relation 
constraints, and represents the primary contribution of this work. In addition, we include 
in this analysis nonzero lepton masses, anticipating future measurements of processes such 
as B — > D ( ~* ) rv T or B — > 7rrz/ r . 

This paper is organized as follows. In Sec. we define the form factors for the physically- 
observed semileptonic decays and present expressions for differential decay widths including 
lepton masses. Section |TJ presents a brief review of the dispersive method for obtaining 
bounds on the form factors. Ingredients of the form factor parameterizations, including ex- 
plicit formulas for the QCD results at one and two loops as well as leading nonperturbative 
effects, appear in Sec. |V|, followed by a tabulation of the form factor weighting functions, 
which are central to the dispersive method. In Sec. |V| we discuss the inclusion of additional, 
previously ignored hadronic states into the dispersion relations, and show how heavy quark 
effective theory may be used to include their effects and thereby tighten form factor con- 



4 



straints. Section [VI] defines and presents the truncation errors that measure the quality of 
our form factor parameterizations. In Sec. |V11| we discuss slope-curvature relations derived 
by Caprini and Neubert [[HJ, point out an invalidating assumption, and examine the form 
of the corrected relations. In Sec. |V111| we present the results of the current analysis in fits 
to the latest experimental data, and in Sec. ITX] we conclude. 



II. FORM FACTORS 

We begin by defining form factors for the semileptonic decays of interest in terms of 
hadronic transition matrix elements. From the field theory point of view, it is most conve- 
nient to define form factors as coefficients of independent Lorentz structures appearing in 
the matrix element. However, the combinations of form factors most easily obtained from 
data are those appearing in a sum of squares in the differential rates, namely, the helicity 
amplitudes. As seen below, the helicity amplitudes are particular linear combinations of the 
original form factors, and thus simply form a different basis for the description of the matrix 
elements. It is these helicity amplitudes we wish to constrain. 

The notation used throughout the paper is as follows: The generic b — > c semileptonic 
decay is denoted by —>■ HJtP, where the hadron iff, has mass M and momentum p, the 
daughter hadron H c has mass m and momentum p', and the charged lepton I has mass 
mi. The momentum transfer t = (p — p') 2 is the invariant mass-squared of the lepton pair 
(or virtual W). The polar vector and axial vector flavor-changing currents are denoted by 
yti = QryU}) anc [ j^n = C7^7 5 6. Finally, it is convenient to define the kinematic invariants 

t± = (M±m) 2 , 

V t V At y J 

where p is the three-momentum of H c in the rest frame of H^. Note that the identity of the 

decaying quark does not enter into the expressions below except through hadron masses, so 

the same expressions apply to such decays as B — > tt£u, D — > K*£u, and so on; only the 

spins of the hadrons are relevant. 
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The inclusion of charged lepton masses brings into the differential widths a new set of 
helicity amplitudes, and thus a new laboratory for studying the strong interaction. The con- 
tributions of such terms is suppressed by a factor of mj, and arises through a virtual W with 
the quantum numbers of a scalar, i.e., a timelike polarization. Angular momentum conserva- 
tion forbids the decay of such a state to a right-handed antineutrino and a left-handed mass- 
less charged lepton, so such decays are necessarily accompanied by a helicity-suppression 
factor mj. Although the detection of lepton mass-dependent effects is presently beyond the 
means of current experiments, the gradual accumulation of statistics may eventually make 
such effects discernible in decays like B — > D^tP or B — > ixtv. 

Finally, it should be noted that the differential widths dT/dt presented below have al- 
ready been integrated in lepton energy, or equivalently over angles of final-state particles. 
This reflects only the current thrust of experiment and does not indicate a limitation of the 
dispersive method described below. Indeed, were the statistics available, it would be inter- 
esting to consider the double differential decay distribution d 2 T/dEidt, for then one could 
probe the parity- violating interference terms between vector and axial vector weak currents. 
In terms of the QCD side of the dispersion relations described below, one would also need 
to compute V-A correlators, in addition to V-V and A-A. 

A. B — ► DIP 

The hadronic matrix element governing the rate of B — > DIP may be described by form 
factors 



(D(p')\V»\B(p))=f + (p + p 



>r+f-(p-p'Y , 



(2.2) 



that enter the differential rate as 



dT _ G 2 F \V cb \ 2 k 
dt 192tt 3 M 3 tl 



\(t- mjf AkH (2t + mj) |/ + | 2 + 3mj\f 



(2.3) 



where 



fo(t) = (M 2 



m 2 )f + (t)+tf_{t). 



(2.4) 
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B. B -> Z>*£i? 



The matrix elements for B — > D*£u depend on four form factors, 



(ir(l/,e)\V*\B(p))=iger*e*JfPj, 

(D*(p\ e)\A»\B(p)) = + (e* • p) [a+(p + pT + a_(p - p')1 , 
that enter the differential rate in the combinations 



where 



dT 
~dt 



G 2 F \V cb \ 2 k 
192tt 3 M 3 t l 



(t-m 2 )' 



(2t + m 2 ) [2t\f\ 2 + \F l \ 2 + 2k 2 t 2 \g\ 2 ] + 3m 2 k 2 t\F 2 \ 2 \ , 



Fi(t) = - 
m 



Flit) = 



m 



2k 2 ta + (t) - -(t- M 2 + m 2 )f(t) 
fit) + (M 2 - m 2 )a + {t) + ta-{t) 



C. A5 — >■ A c £i/ 

This decay can be described by six form factors, defined by 

(A c (p')\V»\A b (p)) = « c (p') [Fir + W + iV"] u b (p), 
(A c (p')|^|A 6 (p)> = ii c (p') [d/ + + G 3 u'1 75« 6 (p), 

with v = p/MA b and v' = p' /M\ c . This gives 

dT _ G 2 F \V cb \ 2 k , 2 x2 
dt ~ 192tt 3 MM V f m n 

x j(t_ -t) (2t + m^) [2t|Fi| 2 + |#y| 2 ] + 3m| (f+ - f) |F„| : 
+ (*+ - t) (2t + mfj [2*|Gi| 2 + |if A | 2 ] + 3m 2 (i_ - i) |G | 



where 



H v (t) = (M + m)F 1+ 1 -(t + -t)Q + ^), 
H A (t) = (M- m)G 1 - i(f_ -t)(^ + ^j, 

F (t) = (M -m)F 1 + ^(t + M 2 -m 2 )F 2 -^-(t-M 2 + m 2 )F 3 , 

G (t) = (M + m)G 1 - -L (t + M 2 - m 2 ) G 2 + ±- (t - M 2 + m 2 ) G 3 . (2.10) 



III. REVIEW OF THE DISPERSIVE APPROACH 

Constraints on a generic Hb — > H c £i> form factor F(t) are obtained by noting that the 
amplitude for production of HbH c from a virtual W boson is determined by the analytic 
continuation of F(t) from the semileptonic region of momentum-transfer m 2 < t < i_ to 
the pair-production region t + < t. The idea of the dispersion relation is to constrain F(t) in 
the pair-production region using perturbative QCD, then use analyticity to translate that 
constraint into one valid in the semileptonic region ||19|| . A detailed derivation can be found 
in [f?|,|8|,pT| ; here we merely outline the essential elements. 

In QCD, the two-point function]] of a flavor- changing current J = V, A, or V — A 

Hf (g) = \W - q'gn^W) + ^fn L j(q 2 ) = i [ d 4 xe^(0\TJ»(x)J^(0)\0), (3-1) 
q q J 

is rendered finite by making one or two subtractions, leading to the dispersion relations 

Tt 2N dU L j 1 r°° , lmU L r(t) 

T/ 2x 1 d 2 U T r 1 f°° , lmU T r(t) , s 



2d(q 2 ) 2 irJo (t-q 2 ) 3 ' 
The functions x(<? 2 ) may be computed reliably in perturbative QCD for values of q 2 far from 
the kinematic region where the current J can create resonances: specifically, 
{m b + m c )Aq C D <C (ra& + m c ) 2 — q 2 . For b — > c, or u, q 2 = satisfies this condition. 



1 This definition differs slightly from that used in [7|HTl||, and serves to separate II ^ into manifestly 
spin-0 and spin-1 pieces. Then the functions x 1 \ X L defined in ( p.2[ ) coincide with x> X defined 
in the previous works. 
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Inserting a complete set of states X into the two-point function relates the Uj to the 
production rate of hadrons from a virtual W, 

IniII^= W(2Tr)W(q-p x )\(0\J\X)\\ (3.3) 
z x 

where the sum is over all hadronic states X with the same quantum numbers as the current 
J, weighted by phase space. Then, from the dispersion relations (|3.2|) , the perturbatively 
evaluated x(g 2 ) is equal to the integrated production rate of W* — > X, weighted with a 
smooth function of momentum transfer squared, t. Since the sum is semipositive definite, 
one may restrict attention to a subset of hadronic states to obtain a strict inequality. In the 
case of interest, we focus on X being two-particle states of the form Hi>H c . This places an 
upper bound on the form factor F(t) in the pair-production region that takes the form 

nx 1 (<r) Jt+ (t - q 2 y 

from the Im IT T dispersion relations in Eq. fl3.2| ). Here W(t) is a computable function of t 
that depends on the particular form factor under consideration. A similar result holds for 

n L . 

Using analyticity to turn ( J3.4j ) into a constraint in the semileptonic region requires that 
the integrand is analytic below the pair-production threshold t <t + . To do this, we introduce 
a function 

z(t; t s ) = *' ~1_= 5 (3.5) 

that is real for t s < t+, zero at t = t s , and a pure phase for t > t + . All the poles in the 
integrand of Eq. Q3.4|) can be removed by multiplying by various powers of z(t] t s ), provided 
the positions t s of the sub-threshold poles in F(t) are known. Each pole has a distinct value 
of t s , and the product z(t;t s i)z(t;t S 2) ■ ■ ■ serves to remove all of them. Such poles arise as 
the contribution of B c resonances to the form factor F(t), as well as singularities in the 
kinematic part of the integrand. After determining these positions phenomenologically, the 
upper bound on F(t) becomes 



/ dt 


dz(t; t ) 


dt 


)t+ 



\<j>(t;t )P(t)F{t)\ 2 <l, 



(3.6) 



where the weight function <p(t; to) (known as an outer function in complex analysis) is given 
by 

W(t) 



<f>(t;t ) = P(t) 



\dz(t;t )/dt\x T (q 2 )(t-q*r 



(3.7) 



The factor P(t) is a product of z(t;t s )'s and y z(t; t s )'s, with t s chosen to remove the sub- 
threshold singularities and cuts in the kinematic part of the integrand, while the Blaschke 
factor P(t) is a product of ,z(t;t p )'s with t p chosen to be the positions of sub-threshold 
poles in F(t). The functions </>(t;to) and P(t) also depend on q 2 , which we leave implicit for 



notational simplicity, while to is a free parameter to be discussed in Sec. |TV. 

The quantity <p(t;to)P(t)F(t) may be expanded in a set of orthonormal functions that 
are proportional to powers of z(t;t ) [see (|4.21|) below]. The function z(t;t ) has a physical 
interpretation as a natural scale for the variation of F(t) in the semileptonic region ||11|| , 
and will play an important role throughout this paper. We exhibit all relevant expressions 
in terms of both variables, z and t, in Sec. |IV|. The result of the expansion in z(t; to) is an 
expression for F(t) valid even in the semileptonic region, 

1 



F(t) 



P(t)<j>{t; h) ^ 



(3.8) 



where, as a result of Eq. ( |3.6| ), the coefficients a n are unknown constants obeying 

oo 

E a l < 1 • (3-9) 

n=0 

For the h — > c processes that are the main subject of this paper, z(t; to) is no larger than 0.07 
for any physical momentum transfer m 2 < t < t_, and can be made substantially smaller by 
a judicious choice of to, so the expansion can be truncated after the first two or three terms. 



IV. PARAMETERIZATION INGREDIENTS 



Generating a parameterization like fl3.8| ) for a particular form factor requires three ingre- 
dients: One needs the perturbative evaluation of \ derived from the two-point function for a 
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current J, including the Wilson coefficients of phenomenologically determined condensates. 
In addition, the functional form of the weighting function <fi must be computed. Finally, the 
masses of sub-threshold resonances with the same quantum numbers as J must be extracted 
from experiment or potential models. The function cf) depends on the form factor under 
consideration, while \ an d P depend only on the current J. 



A. QCD Evaluation of x 

In the previous section we observed that it suffices to take q 2 = in the computation of 
Xj' L for currents containing a heavy quark. This is convenient since then the perturbative 
expressions become particularly simple. Corrections to the perturbative result may be in- 
cluded by expressing the two-point function as an operator product expansion (OPE) and 
including the leading nonperturbative vacuum condensates such as (G 2 ) and (qq); the total 
X is the sum of the perturbative and condensate terms, 

X = Xpert + Xcond- (4.1) 

The full perturbative expressions to two loops may be obtained through a lengthy but 
straightforward 2 manipulation of results existing in the literature [p0| , ^| . At q 2 = 0, x per 



t 



is only a function of the ratio of quark pole masses u = m c jm^^ and for a vector current 
J = V is given by 



327T 2 (1 - W 2 ) 



[1 - u 2 )(3 + Au - 21m 2 + 40u 3 - 21m 4 + 4m 5 + 3m 6 



2 We have corrected a number of typographical errors in the perturbative results of Ref. [20| and 



the Wilson Coefficient results of [ 21 ] to ensure compliance with various consistency conditions. We 
also made several comparisons between references to verify their agreement, once these corrections 
are made. 
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+12m 3 (2 -3m + 2m 2 )1ii(m 2 



a, 



576tt 3 (1 -m 2 ) 6 

"(1 - m 2 ) 2 (75 + 360u - 1031m 2 + 1776u 3 - 1031m 4 + 360m 5 + 75m 6 ) 
+4m(1 - m 2 )(18 - 99m + 732u 2 - 1010u 3 + 732m 4 - 99u 5 + 18u 6 ) 1ii(m 2 ) 
+4m 3 (108 - 324m + 648m 2 - 456m 3 + 132u 4 + 59m 5 - 12m 6 - 9u 7 ) ln 2 (u 2 ) 
+8(1 - m 2 ) 3 (9 + 12m - 32m 2 + 12u 3 + 9u 4 )Li 2 (l - u 2 ) 



(4.2) 



Xpert ( U ) 



7T 2 (1 -U 2 ) 3 



1 - M 2 ) (1 + M + M 2 ) (1 - 4M + M 2 ) - 6M 3 lll(M 2 ^ 

1 - m 2 ) 2 (1 - 36m - 22u 2 - 36m 3 + m 4 ) 



48tt 3 (1 -u 2 ) 4 

-2m(1 - u 2 )(9 + 4m + 66m 2 + 4m 3 + 9m 4 ) 1ii(m 2 ) 

-4m 3 (9 + 18u 2 - 2u 3 - 3m 4 + m 5 ) 1ii 2 (m 2 ) + 8(1 - u 2 ) 3 (l - 3m + m 2 ) 

xLi 2 (l -m 2 )1, 



where the dilogarithm is defined by 

U 2 {z) = 



dz 



, ln(l - z') 



10 Z' 

Expressions for an axial current J = A are obtained from 



!3|) by replacing u 



(4.3) 



(4.4) 



-u. 



It has been pointed out [23| that non-analytic quark mass dependence, such as in the 
perturbative results presented above, indicates the inclusion of some infrared effects into 
the Wilson Coefficients (WC's), in conflict with the usual interpretation of the OPE as a 
separation into short- and long-distance effects. A formal analysis reshuffles the WC's in 
front of each nonperturbative condensate. However, since our analysis requires only the 
numerical sum of such effects, the total result should be the same in either form. 

The leading nonperturbative corrections are supplied by the condensates of dimension 
less than five, namely the gluon and quark condensates. For a heavy h quark decaying into a 
quark q of arbitrary mass through a vector current, the contributions from the condensates, 
derived using Refs. ||2~1~1 , |2"4| and evaluated at q 2 = 0, read 



12 



xLd(w) 



. (2-3m + 2m 2 ) 
■(??>■ 



2m\(\ - a 



2\5 



+ 



X 



S G 2 



,2V 



v 24mf(l - u , 

1 - m 2 ) (2 - 104m + 148m 2 - 270m 3 + 145m 4 - 104m 5 + 5m 6 - 2m 7 ) 

-12w 1ii(m 2 ) (3 - 5m + 17m 2 - 15m 3 + 17m 4 - 5m 5 + 3m 6 )] j , 



(4.5) 



Xcond( M ) 



+ (gg) m 3 n -m 2 ) 3 



a. 



+ (^G 2 



12m 4 (l -m 2 ) 5 



(l - u 2 ) (l - 21m + 10m 2 - 20m 3 + 
-3m1ii(m 2 ) (3 - 2m + 8m 2 - 2m 3 + 3m 



M 4 -M 5 



These expressions superficially appear to diverge in the limit u 
limit the {m q qq) condensate obeys the relation 



(4.6) 



1. However, in this 



(4.7) 



and Eqs. (|4.5| ), (|4.6| ) simplify to 



xLd(w) 



a 



s ^2 



vr / 24(1 -M 2 ) 7 m^M 
x j (l - m 2 ) (2 - 5m + 104m 2 - 145m 3 + 268m 4 - 145m 5 + 104m 6 - 5m 7 + 2m 8 ) 

+12u 2 ln(M 2 ) (3 - 5n + 17m 2 - 15u 3 + 17m 4 - 5m 5 + 3m 6 ) j , (4.8) 



Xcond( U ) 



■G 2 ) 1 . — jfl - u 2 ) ( 1 - m + 20u 2 - 10m 3 + 20m 4 - m 5 + m 6 ) 

/ 12(l-M 2 ) 5 m 4 M\ V A > 

+3m 2 ln(M 2 ) (3 - 2m + 8m 2 - 2m 3 + 3m 4 ) j , (4.9) 



The expressions ( |4.8|) , (|4.9|) are applicable to b — > c decays, while (|4.5|) , ( |4.6|) are best 
suited to b — ► u or c — ► s transitions. Expressions for an axial current are obtained by 
replacing u — > —u. The limiting values for u — » are 
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mix 1 (0) 

x L (o) 



3 a s . , 2 \ 1 , . I 



+ -^-(25 + 4^)-^^) 



32vr 2 192vr 3 ' my, 
87r 2 1447T 3 mi 



12m£ \ 7T 



+ 



a 



while for u — > 1, 



™ 2 x T (-i) 



1 41a, 
+ 



20vr 2 162vr 3 210m£ 
1 689a, 1 



12m£ \ vr 



— G 2 

7T 



G 



40vr 2 6480vr 3 140m 4 \ it 



a 



G 1 



1 7a s 1 



4tt 2 12vr 3 60m£ \ vr 



a 



s n2 



For 6 — > c, using pole mass values such that k = 0.33, we have 



m^ T (+0.33) = 9.659 ■ 1(T 3 
m^ T (-0.33) = 5.709 ■ 10~ 3 
X L (+0.33) = 3.713- 10~ 3 
X L (-0.33) = 2.162- 10~ 3 



1 + 1.42a, - 4.8 • 10 
1 + 1.32a s - 6.8 ■ 10 
1 + 1.37a, - 5.3 ■ 10 
1 + 0.64a s + 2.2 ■ 10 



-4 



4.9GeV\ 4 / (a s G 2 /n) ' 
) \ 0.02 GeV 4 



m b 



4.9GeV\ 4 / (a s G 2 /n) 



m b 



-) 



0.02 GeV 



-4 



4.9GeV\ 4 / (a s G 2 /n) 



m b 



0.02 GeV 4 



. t ,4.9GeV\ 4 / (a 8 G 2 /ir) 



0.02 GeV 



(4.10) 



(4.11) 



(4.12) 



where the central value < ^-G 2 >= 0.02 GeV 4 is taken from Ref. |26| , and the pole mass 
value m b = 4.9 GeV is from f27|. We also use a s {m b ) = 0.22 [p7H , and since the coefficient 
of the gluon condensate is tiny, we ignore it in our numerical analysis for b — > c. 



B. Weighting functions <fi 

To obtain the general form of the weighting functions (f) defined in Eq. ( |3.6|) , first observe 
that the quantities U L , n T are, respectively, the /i = v = and fi = v = 1, 2, or 3 components 
of H 7 ^ evaluated in the center of mass frame, = (^/t, 0) . Then the generic expression for 
the contribution of a particular form factor F(t) to the polarization tensor may be denoted 
by 
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ImIT T > ^-{t - f + )t(t - «_)* t- c \F(t)\ 2 6(t - t+), 

lmU L > —(t - £ + )f (t - r( c+1 )|Fmi 2 0(f - t + ), (4.13) 

where if, a, b, and c are integers determined by the form factor F(t), and rij is an isospin 
Clebsch-Gordan factor, which is 2, 3/2, and 1 for B — > B — > 7r, and A;, — > A c transi- 
tions, respectively. Also, let x — X T or X L denote the generic QCD function appropriate to 
the quantum numbers of the form factor] 3 ] The expressions for the weighting functions are 
readily derived from Eq. Q4.13 ), and are given by 



/ 711 


f t+-t\ 


1 Kn X 


\t+-to) 



<Pi(t; to) — J j7~~ f ~j~~r ) (v^+V^-M (*+-*) ; 



f / r. r r—\ -( c + 3 ) 



x{y/t + -t + yjt + -t_) (yjt + -t+y/t+) . (4.14) 

The values of the parameters K, a, b, c for each form factor, as well as the relevant \i are given 
in Table 1 for B — > transitions, and in Table 2 for A& — > A c . Although ( (4.14| ) assumes 
q 2 = 0, it is easy to generalize to arbitrary Q 2 = —q 2 : Simply evaluate the perturbative 
functions x{Q 2 ) ari the given value, and multiply (f> by 

. / v d 

u-t + . 

(4.15) 



/ t+-t + yjt+ + Q 2 
with d = 3 if the form factor involves \ T and d = 2 if it involves x L - 



3 In terms of previous notations, one finds for the meson form factors considered in Refs. 
K = 2 p k~ 1 /irxM 2 , a = b = p, c = s + p — 3, whereas in Ref. Jll| one finds K = 3-2 s ,a = b = w, 
and c = p. For the baryon form factors considered in M, the relation is given by by a = 2p + 1, 
6 = 3 — 2p, c = s + 1, and if = 2k~ 1 /-kxM 2 or K = 2k~ 1 / ttx > depending upon whether the form 
factor appears with x or x L . 
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Ft 


K 


X 


a 


b 


c 


u 


48 


x T 

/V 


'+u) 


3 


3 


2 


fo 


16 


L 

A. 


'+u) 


1 


1 


1 


f 


24 


x T ( 


—u) 


1 


1 


1 




48 


x T ( 


—u) 


1 


1 


2 


9 


96 


x T 


» 


3 


3 


1 


T 2 


64 


x L ( 


—u) 


3 


3 


1 



Table 1. Factors entering Eq. ({1.14 ) or ( 4.23 ) for the meson form factors Fi in B — > 



Fi 


K 


X 


a 


b 


c 


Fo 


8 


x L 


[+u) 


3 


1 


1 


F x 


12 


T 

X 


[+u) 


1 


3 


1 


H v 


24 


T 

X 


[+u) 


1 


3 


2 


Gq 


8 


x L 


[—u) 


1 


3 


1 


Gi 


12 


T 

X 


—u) 


3 


1 


1 


H A 


24 


T 

X 


—u) 


3 


1 


2 



Table 2. Factors entering Eq. Q4.14\ ) or ( 4.23 ) for the baryon form factors Fi 



While the momentum-transfer variable t = (p — p') 2 is useful for heavy-to-light decays 
and has an obvious physical meaning, it is often more convenient when dealing with heavy- 
to-heavy transitions such as b — >• c to use a kinematic variable that helps disentangle long- 
distance physics from the heavy quark scale. One such variable is 

, M 2 + m 2 -t 



w = v ■ V 



(4.16) 



2Mm 

In the b rest frame, w depends only on the energy transfer to the light degrees of freedom, 
in units of Aqcd- It is due to this property that b — > c form factors are related to each 
other in the heavy quark limit at equal values of w. In the semileptonic region, the variable 
z(t; to) has the same physical property, so form factors will be related at equal values of z 
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as well. We can demonstrate this by expressing z as a function of w without reference to 
heavy meson masses, 

z(w; N) = \ = , (4.17) 

where iV is a free parameter related to t by 

N= f f^, (4.18) 



so that z(w;N) vanishes at w = 2N — 1. With t,to related to w, N by Eqs. ([4.1 6|) and 



( |4.18| ), Eq. (|4.17| ) is simply a rewriting of Eq. ( |3.5|) , z(t;t ) = z(w;N). The advantage in 
using z(w; N) for b —>■ c transitions is that its definition is process- independent. 

For the semileptonic decay H b — > H c £v, the limiting values of z are given by mj <t<t_, 



or 



.y/N+ 1 

and 



^ = /^zl) (4 . 19) 



1 + r )2 _ p _ 2 ^Nr 



(4.20) 



y l+r) 2 - 5 2 + 2VNr 

where r = m/M and 5 = m^/M. 

The dispersion relation (|3.6|) , written now entirely in terms of z, reads 

/ -|0(^P(,)F(^)| 2 <1, (4.21) 

where C is the unit circle in the complex z plane, the Blaschke factor for a pole at z p = 
z {tp] to) (which is real for sub-threshold resonance masses) is 

z(t;t p ) = ^^, (4.22) 
1 — zz p 

for z = z(t; t ) and any t , with P(z) being the product of all such factors, and the weighting 
functions ft4.14|) are given by 
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V Kir x 



x 



-(c+3) 



(4.23) 



(l + r)(l- z) + 2VNr(l + z) 

Finally, evaluation of 4>(z; N) at nonzero values of Q 2 = —q 2 is accomplished by multiplying 
(H)by 

[l + r)(l~z) + 2VW(l + z) , 



r (l + r) 2 + Q 2 /M 2 (l - + 2ViVV(l + 2) 
where d is defined as in Eq. ( [4.151) . 



C. Sub-threshold Resonances and P(t) 



The Blaschke factor P(t) for a form factor describing Hf, — > if c £i/ depends on the masses 
of 5 C resonances below the Hj,H c pair-production threshold. The Blaschke factors are simply 
products of z(t] t a ) with t s evaluated at the invariant mass-squared t s 
resonance with the same spin-parity as the current J, 



t p of each such 



p(t) = II*(*;^)- 



(4.25) 



For b — > c transition from factors, the masses of the relevant -B c -type resonances can be 

We compile in Table 3 the masses 



accurately estimated from potential models [28 



computed in Ref. [29 



Type 


Masses (GeV) 


Vector 


6.337, 6.899, 7.012, 7.280 
7.350, 7.594, 7.646, 7.872, 7.913 


Axial vector 


6.730, 6.736, 7.135, 7.142 
7.470, 7.470, 7.757, 7.757 


Scalar 


6.700, 7.108, 7.470, 7.757 


Pseudoscalar 


6.264, 6.856, 7.244, 7.562, 7.844 



Table 3. Calculated B c pole masses used in this work. 
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For heavy-to-light form factors there is no formal limit in which the light quark becomes 
nonrelativistic, and potential model calculations are less reliable. This is not a problem for 
B — > 7r, where the only sub-threshold resonance, the B*, is experimentally observed, but 
for decays to other light states such as B — > p, the presence and masses of additional sub- 
threshold resonances must be taken from models. Once these uncertainties are accounted for, 
simple parameterizations should be reliable. For example, the model of indicates that the 
form factor / for B — > p has only one narrow sub-threshold pole. Indeed, this pole appears 
to have been observed (mixed with others) by ALEPH, DELPHI, and OPAL [see |27] for 
analysis and references on the "B}(5732)"]. Even accounting for significant uncertainties in 
its mass, this leads to an accurate parameterization using the overall normalization and two 
unknown coefficients. It is important to estimate the uncertainties from model-dependent 
poles on a case by case basis |L4[ . 

This applies as well to sub-threshold branch cuts due to multi-particle states and anoma- 
lous thresholds. A model-dependent analysis 0|| suggests these are negligible for b — > c 
transitions. Qualitatively, this result comes about because cuts are a much less severe form 
of non-analytic behavior than poles. Whether cuts continue to be unimportant for B — > p 
transitions requires a more detailed analysis. 



V. ADDITIONAL STATES 

The effects of higher states in the dispersion relation depend on the flavor of the b — > q 
current under consideration. Henceforth, we specialize to b — > c transitions, for which HQET 
is most useful. 



A. Contributions to the Dispersion Relation 

We have observed that contributions to the original dispersion relation [see Eq. ( |3.3| )1 are 
semipositive definite, so each additional state coupling to the vacuum through the current 
J serves to further saturate the bound supplied by the QCD parton-level calculation. The 
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inclusion of only a single two-particle BD or BD* state in obtaining these bounds is relatively 
weak, since such exclusive states account for only a small portion of the inclusive total. In 
general, the hadronic side includes also B c resonances, a continuum of states like BD7nr, 
and so on. While it is desirable to include as many of these states as possible, it is not 
clear how to include them in a model-independent fashion; the chief exceptions are two- 
particle states related to one another via heavy quark spin symmetry, namely the four 
states B(*'D(*\ While form factors for transitions such as B* — > D^*' are not physically 
accessible through semileptonic decays, their normalization is nonetheless known via HQET, 
allowing an additional strengthening of the dispersive bounds. 

The contributions ( 4.13Q of these states to the dispersion relations are of the form 



ImII(t) >E I 2 » O 5 - 1 ) 

i 

where the sum is over all helicity amplitudes Fi arising from pair production of B^D^*'. 
The weight functions (f> are readily obtained from the kinematic prefactors while the 

relation between the helicity amplitudes and the original form factors may be obtained by 
choosing definite polarizations of the B*, D*, and the virtual W. 
For B* — > D transitions, the form factors are defined by 

(D(p')\V^\B*{p,e))=ige^e a p p p' T 

(D*(p')\A»\B*(p,e)) = /> + (e • p') [a+tf + pY + a_(p' - p)»] . (5.2) 

The functions <fi for B* —>■ D transitions are identical to those for B — > D*, with the simple 
replacement M <-> m. The helicity amplitudes possess factors of M and m, and need not be 
invariant under this exchange. This is true in particular for T\ and JF 2 ; here we find that 
[compare ( |2.7|) 1 



2k 2 ta + {t) - i(t - m 2 + M 2 )f(t) 



f(t) - (M 2 - m 2 )a+(t) + ta-{t) 



(5.3) 
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The transition B* — > D* possesses 10 independent vector current form factors, which we 



define as follows pT | 
1 



rn 



(D*{e D ,v')\V»\B*{e B ,v)) = U(e B ■ e* D )v'» + f 5 (e B ■ e* D W 



+UW ■ e B )(v ■ e* D )v» + f 7 (v' ■ e B )(v ■ e* D )v'» 
+fs(v-e* D )e B + f 9 (v'-e B )e^, 



VMm 



(D*(e D ,v')\A»\B*(e B ,v)) = +if w e» aM e B e*gv^ + z/ne^e^V 

+if 12 [W ■ e^e^e^VV + e af3yS e B e*£v"tv 5 v'» 
+*/i3 \{v ■ e D )e^ af3y e* B a v'^ - e^e^V^V 



(5.4) 



The combinations of these form factors appearing as helicity amplitudes may be denoted 
/ 9 



V ++ 

So+ 

A + 
A+o 



h 
1 

\/2Mm 



(MU + mf 5 ) 



— J— T {Mm (t - M 2 - m 2 ) (M/ 4 + mfe) - 2k 2 t (mf 6 + Mf 7 ) 



+Mm 
1 



(2Mm) 
+2k 2 t 
1 

VIMm 



m (t + M 2 - m 2 ) / 8 + M (t - M 2 + m 2 ) / 9 ] }, 
f {Mm (t — M 



2 2 

— m 



M (t — M 2 + m 2 ) f 4 -m(t + M 2 - m 2 ) / 5 
(t + M 2 - m 2 ) f 6 -M(t-M 2 + m 2 ) / 7 ] - Ak 2 tMm (m/ 8 - M/$ 
M(t- M 2 + m 2 ) f 4 -m(t + M 2 - m 2 ) / E 



~~77~ nr {Mm [(t — M 2 — m 2 ) / 10 - 2Mm/nl - 2A; 2 t/ 12 } 
2(Mm)2 >■ lv / J 

— r (M 2 m ft - M 2 + m 2 ) f 10 - Mm 2 ft + M 2 - m 2 ) f 

2(Mm)5 1 v y v y 

+2A; 2 t (M/ 12 -m/ 13 )}, 

f {Mm [2Mmf 10 - (t - M 2 - m 2 ) / n ] + 2A; 2 t/ 13 } 



n 



1 



2(Mm)2 
1 



2(Mm) 



f {2Mm (Mf w +mf n ) + [m (t- M 2 +m 2 ) f 12 + m (t + M 2 - m 2 ) / 13 ] } 



(5.5) 
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The labels V, S, A, P reflect the spin-parity (vector, scalar, axial vector, or pseudoscalar) of 
the virtual W, while the subscripts denote the helicities of the W* and D* in the decay of 
the B*. 

Only amplitudes of a fixed spin-parity enter each dispersion relation. For example, the 
four V helicity amplitudes enter the vector current dispersion relation for n T . Above all the 
appropriate pair-production thresholds, the contribution to IulIT^ from BD, B*D, BD*, 
B*D*, and A;,A C states is 

k 3 



2(\f + ^+\V 0+ \ 2 + \V 00 n+t(\g\' + \g\' + \K 



\V + o\ 2 ) 



+ 



k(t - t_) 

127rf 3 / 2 



2t|F 1 | 2 + \H V \ 2 



(5.6) 



where we have included an isospin factor rij = 2 for the mesons, and k is defined in Eq. ( |2.1| ). 
Using Eq. ( |4.13| ), it is straightforward to compute the weighting functions cf) for all the B*D* 
form factors. For each such Fj, one obtains a parameterization of the form Eq. ( |3.8| ), with 
the unknown expansion coefficients denoted by 6j n , 

1 



FAz] 



W n=0 



(5.7) 



The <p parameters [see ( |4.14 ) or ( 4.23 )1 a i b, c, K and the relevant \ for B* — > D* are given 
in Table 4. 
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8 


Y L 




1 


1 


1 
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24 


T 

X 


:-«) 


1 


1 


1 


A 0+ 


24 


T 

X 


:-«) 


1 


1 


2 
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24 


T 

X 


:-«) 


1 


1 


1 


Po+ 


32 


x L 


:-«) 


3 


3 


1 



Table 4. Factors entering Eq. ($.14 ) or ( ]4.23| ) for the meson form factors Fi in B* — > D 



Substituting this expansion into the dispersion relation ( 4.21|) gives 



H oo 
i=0 n=0 



(5.8) 



Included in the sum are all helicity amplitudes i = 0, • • • , H for processes with the right 
quantum numbers to couple to the current J. It is clear that the constraint on a particular 
helicity amplitude Fi can be strengthened if it is possible to relate the various bi n , i.e., if one 
can relate the form factors. This is accomplished with the help of heavy quark symmetry. 



B. Form factors in the Heavy Quark Limit 

In general, form factors are not related by heavy quark symmetry throughout the pair 
production region |j32| , but may be related in the semileptonic region. This fact has been 



exploited to improve constraints on the B — » B elastic form factor |17| and the dominant 
B — > 7i£u form factor |Tl|]. The situation in the present case is conceptually analogous, 



although algebraically more cumbersome. 
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In the semileptonic region, the 20 form factors of — > £)M reduce to only one in the 
heavy quark limit, the universal Isgur-Wise function £(«;), with £(1) = 1. Likewise, the six 
baryon form factors in A& — > A c reduce to another universal function |33j| , which we may 
denote C( w ); with ((1) = 1. Recalling that r = m/M, the relation of the helicity amplitudes 
to £ and £ in the heavy quark limit are given by 

/+ - -f - -ft - -^v 0+ - v 00 - -p 0+ - -^t, 

f = jr i = jr i = -^2S 00 = S 0+ = -A 0+ = M 2 v^(l - r)(l + w)£, 
9 = -9 = V ++ = V +0 = ^, 

f = -f = ~A ++ = A +0 = M^(l + w)Z, (5.9) 
for mesons, and by 

F = H A = M(l-r)(, 
Fx = G x = C, 

H v = G = M(l + r)( (5.10) 



for baryons. The only helicity amplitudes in this list protected by Luke's theorem PJ5^ 
from 1/M corrections at w = 1 are /, /, and G*i. 

In the strict heavy quark limit, our dispersion relation constraints become useless because 
the form factors in Eqs. ( |5.9|) and (|5.10|) develop an essential singularity due to an infinite 
number of poles just below threshold ||35|| . The description Eq. ( |3.8| ) then contains no 
information, because the Blaschke factor P(t) goes to zero in the semileptonic region. 

For finite masses our parameterizations are well-behaved, and heavy quark relations are 
valid up to 1/m corrections. In this case, some of the form factors in Eqs. ( |5.9| ) or (|5.10| ) will 



have more constraining parameterizations than others because their Blaschke factors, which 
reflect the number and positions of sub-threshold resonances, will be larger. If one views 
the universality of the Isgur-Wise function as arising from the dominance of the essential 
singularity in each of the form factors, one might expect 1/m corrections to be larger for form 
factors with fewer sub-threshold poles. Of course, it is always possible that the residues of 
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the various poles could conspire to keep 1/m corrections small, since Blaschke factors alone 
permit |36[ any residues consistent with the dispersion relation bounds. 

We use the heavy quark relations ( |5.9| ) and ( |5.10| ) to tighten our bounds on parameter- 



ization coefficients. When we use these bounds to quote smaller errors on our parameter- 
izations, we allow for substantial deviations due to 1/m effects, thereby minimizing errors 
induced by assuming full heavy quark symmetry. However, when we use these bounds to test 
heavy quark symmetry by constraining the slope or curvature of the Isgur-Wise function, 
one should bear in mind the possibility that heavy quark violations could be larger in form 
factors with very few sub-threshold poles (such as fo, Sqo, or 5*o+) than in those that are 
typically measured experimentally. 

C. Bounding Parameterization Coefficients 

We now use the heavy quark relations of the higher resonance helicity amplitudes to 
improve the constraints. For concreteness, consider the form factor g(w). Near zero recoil, 
w = 1, heavy quark symmetry relates it to six other form factors appearing in the dispersion 
relation, 



g (w) = -g{w) = V ++ = V-i 



o 



2 ~f+( w ) = TTl = ~T7~T~ — ^o+- ( 5 - n ) 



M + m M + m M + m 

While exact heavy quark spin symmetry implies that the functional dependence of these 

form factors is the same for all values of w corresponding to semileptonic decay, for physical 

masses we demand only that the normalization and first and second derivatives of these form 

factors are roughly equal at w = 1. This suffices to provide lower bounds, assuming full 

heavy quark symmetry, on the contributions of the form factors in Eq. ( [5.11| ) to the sum in 



Eq. ( |5.8| ). Once computed, it is straightforward to weaken these lower bounds by including 
factors indicating the violation of heavy quark symmetry. 

Since each form factor F{ in ( p.ll| ) has a parameterization of the form ( p.7[ ), the coeffi- 
cients a n in the expansion of g can be related to the coefficients bi n in the expansion of Fi 
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by 

E^^rnE 6 *.^ ( 5 - 12 ) 

where ~ means only the normalization and first two derivatives are equal at w — 1. The 
functions 

Q = ^#^ (5.13) 

are given by ratios of Blaschke factors, weighting functions, and symmetry factors 
chosen so that 

C{z) = EitfFiiz) (5.14) 

in the heavy quark limit. The kinematic factors Hj(w) appear in Eqs. (|5.9| ) and ( |5.10| ). By 
choosing the same value of N for each form factor, we ensure that the kinematic variable 
z(w;N) is process-independent, most of the z dependence in <f> g (z) /<f>i(z) cancels out, and 
the Ci(z) become quite simple. 

Numerically, the values of N that optimize our constraints correspond to iV = 1 + 2e 
with e ~ 0.05, and ignoring terms of order e 2 in Eq. (|5.12 ) is a good approximation for 



parameter values a n that saturate their bounds. Evaluating ( 5.12 ) and its first and second 
derivatives at w — 1 (z ~ — e/2), we find for each i, 

b = Ca 

b\ = C'a + Ca\ 

b 2 = \c"a + Cat + Ca 2 + ^e(b 3 - b* QS ) , (5.15) 

where C, C = dC/dz and C" = d 2 C/dz 2 are evaluated at z = (notw = 1), and is the 
value 63 would have if the third derivative of ( |5.12| ) yielded a valid relation. Departures from 
the heavy quark symmetry limit tend to increase as we take higher derivatives of ( |5.12| ), so 



one might expect substantial corrections to 63 = 63 However, corrections to this relation 
are multiplied by e, so we may justifiably ignore the factor 3e(&3 — 6 3 iCJS )/2 in b 2 . 
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Substituting Eq. ( |5.15| ) into 

oo H oo 

E^+EEC<i (5.i6) 

n=0 i=l n=0 

gives more stringent bounds on the coefficients a n . For the form factor g, they yield —0.37 < 
ai < 0.40 and —0.49 < a 2 < 0.47, compared to |ai|, \a 2 \ < 1 from Eq. (|3.9|). Bounds for the 
other form factor coefficients can be obtained in a straightforward fashion by singling out a 
different set of b in in (|5.12j ). 



The bounds on a\ are useful as tests of heavy quark symmetry. For example, parameter- 
izing /o and constraining —0.62 < a\ < 0.58 restricts the slope of the Isgur-Wise function, up 
to l/m c effects, to —0.3 < p 2 = — (d£ / dw)\ w= i < 1.8. Bounds on a 2 are useful for decreasing 
truncation errors of our one-coefficient parameterizations, as described in the next section. 
In this case, however, we do not want to rely on exact heavy quark symmetry, because we 
wish to be as conservative as possible with how much the symmetry improves the trunca- 
tion error, so we allow for explicit violations to the infinite mass limit. Such violations are 
potentially largest for the bounds on a 2 , because they depend on relations involving second 
derivatives. Allowing for the total contribution from the higher-spin states to be as low as 
60% of their infinite-mass value gives bounds on a 2 of —0.58 < a 2 < 0.57 for g, which leads 
to a corresponding decrease in the error induced by truncating the expansion (|3.8|) by 40% 
[see Eq. (|6.4| ) below]. In summary, bounds on a>i are obtained using full, unbroken heavy 
quark symmetry, and may be used to test the accuracy of this symmetry. The bounds on 
ai do not enter the construction of our parameterizations. The derivation of the bounds on 
a 2 allows for substantial violations of heavy quark symmetry. These bounds enter into the 
truncation errors we quote on our parameterizations. 

Since the baryonic and mesonic form factors are not related by heavy quark symmetry, 
Eq. (|5.16| ) only applies in the baryonic case to the pairs Fi,H v and Gi,Ha- Of these, 
only Hy and Ha receive substantial improvements to their truncation errors, which are 
proportional to 03 [see Eq. ( |6.4| )]. While l/m c corrections to relations involving 03 and 63 
could a priori be large, they are known in terms of one constant Aa for the baryons |37 
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and are not particularly large. For example, the relation between Hy and F\ is independent 
of w, to 0(l/m^). Allowing for the contributions to the dispersion relation from i*\ and G\ 
to be as small as 50% of their heavy quark symmetry values gives 1 0.3 1 max = 0.57 for Hy and 
|a3|ma X = 0.29 for H A . 

The bounds on a x ignoring heavy quark violation, and the bounds on a 2 allowing for 
deviations from the heavy quark limit as described above, are given in Tables 5 and 6. 



VI. TRUNCATION ERRORS 

To fit data with our parameterizations, we must truncate the series (p.8|) after a finite 
number Q of unknown coefficients a n . This introduces a truncation error which can be 
minimized by choosing an appropriate value of N [[/J (or equivalently, to)- Rather than use 
ao as a free parameter as in ||, we solve for ao so that the form factor is automatically 
normalized at w = 1 to F(l). If N = 1, then z(l;N) = 0, and this parameterization 
coincides with that in ||. For N > 1, the parameterization ( |3.8| ), including the solution for 
a , may be re-expressed as 

E F (w)F(w) = - I P(1)0(1; N)F(1) + £ a n [z n (w; N) - z n (l; N)} 1 , (6.1) 



P(w)(j)(w;N) 



n=l 



where the product SpF is normalized to coincide with the Isgur-Wise function £(u>) [or 
C(w) for baryons] in the heavy quark limit. The full form factor is of course given by the 
Q — » 00 sum, while approximations F Rt are obtained by truncating at finite Q; then the 
fit coefficients ai, . . . , ag can be chosen so that the difference between the parameterization 
F fit and the actual form factor F is given by 

E F (w) (F(w) - F & \w)) = f f ffi £ a n (z n - z 
V ' P{w)<t>{w) n ^ +1 



mm / 



a Q+1 (^ +1 -^n), (6-2) 



P(w)(j)(w) 

where we have ignored numerically unimportant higher order terms. Using the Schwarz 
inequality, the boundedness condition J2 n a n ^ 1> an d the geometric series sum, these higher 
order terms can be shown to be smaller than 
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^ f a < z r. 



< 



< 



P(w)<f)(w) 
E F (w) 



\ 



J2 K 

n=Q+2 



; £ 1 

n=Q+2 



Z' L - Z„ 



— z, 



Q+2 



+ 0(z 



(6.3) 



P(w)(j)(w) 

which, since \z\ < 0.04 for the cases at hand, can be justifiably ignored. Thus, to good 
accuracy, the maximal possible truncation error T(SpF) for any physically allowed z is 



T(EpF) = max 



(w)(z 



Q+l 



|max; 



(6.4) 



P(w)(f)(w) 

where Q is the highest power of z used in the parameterization fit and |aQ + i| max is the 
maximal allowed value of The inclusion of higher states leads to tighter bounds on 

|flQ+i|max and thus smaller truncation errors. 

The truncation error vanishes at w = 1, where the normalization is fixed. For optimal 



values of N, z min « 



so for one-coefficient parameterizations the truncation error 



drops sharply for some w near w max . This leads to significantly smaller truncation errors 
than in previous works J7]||. For two-coefficient parameterizations, the factor (z 3 — z^ in ) 
adds constructively near w max ; for a given bound on 1 0,3 1 , this leads to a larger truncation 
error than in previous work 0. This is unavoidable as long as only a\ and 02 are fit 
parameters, i.e., as long as ao is chosen to enforce the normalization at zero recoil, so that 
the truncation error vanishes at z — z m \ n rather than z = 0. 

Our truncation errors for one- and two-coefficient parameterizations of the various form 
factors in B — > D£u, B — ► D*£u, and A& — > A c £u are shown in Tables 5 and 6. The 
value of the free parameter N has been optimized for each form factor and number of 
fit parameters to produce the smallest truncation errors. The bounds on a\ come from 
Eq. (|5.16|) and the heavy quark symmetry relation Eq. (|5.15|) , and may be used as tests 
of heavy quark symmetry. The bounds on a 2 enter our truncation errors, so allowance for 
heavy quark symmetry violation has been made as described in the previous section. The 
truncation error is expressed as a percentage of S^(1)F(1), which equals unity in the heavy 
quark limit. The truncation errors in Table 5 are typically less than half those of previous 
parameterizations, while the errors in Table 6 are either better or worse by nearly a factor 
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of two, depending on the form factor (we have corrected an oversight in this reference, 
which used incorrect Blaschke factors for F and Go). 

Note the especially small size of the truncation error for the form factors fo and JF 2 ; if 
we consider all form factors related by heavy quark symmetry, then the champion in this 
respect is Sq + , with a truncation error of only 0.56%. However, as we discuss in the next 
section, actually fitting data to such form factors introduces much larger 1/m uncertainties. 



F 


-^optimal 


Combined bounds from B^D^ 


T(E F F) 


u 


1.108 


-0.23 < ai < 0.20 


-0.55 <a 2 < 0.58 


2.6% 


fo 


1.109 


-0.62 < ai < 0.58 


-0.78 <a 2 < 0.85 


0.7% 


f 


1.093 


-0.37 < ax < 0.39 


-0.58 <a 2 < 0.56 


1.8% 


Fx 


1.093 


-0.06 < ai < 0.06 


-0.11 < a 2 < 0.10 


2.1% 


9 


1.093 


-0.37 < ai < 0.40 


-0.58 <a 2 < 0.57 


1.2% 


T 2 


1.093 


-0.41 < ai < 0.45 


-0.59 <a 2 < 0.57 


0.6% 




1.081 


-1 < ai < 1 


-1 < a 2 < 1 


6.3% 


Go 


1.080 


-1 < a x < 1 


-1 < a 2 < 1 


9.4% 



Table 5. One-coefficient parameterizations using optimized N. Bounds on ai ignore heavy 
quark violation. Bounds on a 2 allow for violation as described in Sec. 0. The truncation 
error is relative to the normalization of SpF at w = 1. 

For baryon form factors, the large number of sub-threshold poles typically ensures that 
at least two parameters are required. Even with two parameters and the spin-symmetry 
improvements, the truncation errors are significant. Using three parameters reduces the 
truncation errors to negligible levels. 
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F 


-^optimal 


Combined bounds from B^*'D^*' 


T(~ F F) 


Fi 


1.104 


-0.90 < ai < 0.90 


-0.97 <a 2 < 0.97 


7.0% 


H v 


1.104 


-0.31 < a x < 0.31 


-0.53 <a 2 < 0.53 


12% 


G x 


1.104 


-0.98 < ai < 0.98 


-0.99 <a 2 < 0.99 


9.1% 


H A 


1.104 


-0.15 < ai < 0.15 


-0.26 <a 2 < 0.26 


18% 


F 


1.104 


-1 < Oi < 1 


-1 < a 2 < 1 


0.35% 


Go 


1.104 


-1 < Oi < 1 


-1 < 2 < 1 


0.53% 



Table 6. Two-coefficient parameterizations using optimized N. Bounds on ai ignore heavy 
quark violation. Bounds on a 2 allow for violation as described in Sec. [V]. The truncation 
error is relative to the normalization of "BpF at w = 1. 

A way to circumvent the relatively large truncation errors on most of the baryon form 
factors is revealed by an interesting feature of Table 5: The A;, — > A c £u form factor Fq is well- 
described (to about 6%) by a one-coefficient parameterization. The contribution from Fq to 
the decay rate is suppressed by the lepton mass, so it is difficult to observe. However, F is 
related by heavy quark symmetry to all the other baryon form factors [see (|5.10|) 1- Unlike 



in meson decays, the relation including 1 / m c effects is known in terms of one constant |37 
Aa ~ Ma 6 — nib. Thus, the differential decay rate can be described in terms of two constants, 
ai and Aa- Using the value of \V c b\ obtained from B — > D(*Hi> and known hadron masses 
should then allow a determination of A w Mb — nib- This quantity is important because it 
enters both exclusive and inclusive semileptonic decay distributions as a l/m c correction. 

In addition to errors incurred by truncating the expansion (|6.1|) , there are a number 
of uncertainties arising from various approximations we have made. We enumerated in || 
a list of uncertainties which must be estimated for a reliable determination of the quality 
of the dispersive bounds, and found that we could allow for their effects by increasing 
truncation errors by 40%. Since then, we have greatly reduced or eliminated many of 
these uncertainties. In 0] we saw that the inclusion of the additional parameter N (or to) 
permits a dramatic reduction of the truncation error, and the more realistic definition of 
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the truncation error used here reduces the uncertainties of one-parameter fits even more. 
The simultaneous inclusion of all B^*'D^ states on the hadronic side of the dispersion 
relation serves to help saturate the bound from the partonic side, and the explicit inclusion 
of two-loop perturbative and leading nonperturbative effects eliminates them as a source of 
uncertainty. The uncertainty in B c pole positions is most significant for poles near threshold, 
and leads to larger truncation errors only if the pole masses have been overestimated. As 
pointed out in 0, branch cuts from multi-particle states below threshold can be ignored 
if they violate isospin. The only uncertainty from || that has not been reduced is due to 
the choice of pole quark masses in the perturbative calculation. Combining the remaining 
uncertainties as in ||, we find that their effects may be allowed for by increasing B — > D, D* 
truncation errors by 20% (e.g., the conservative truncation error for / would be 2.2%) and 
Afc — > A c truncation errors by 30%. In nearly all cases, this small increase makes no practical 
difference. 



VII. SLOPE AND CURVATURE RELATIONS 

A set of interesting relations between the slopes and curvatures of B — > Div and 
B — > D*iv form factors has been derived by Caprini and Neubert [[TJJ. Here we examine 
these relations in the context of the parameterization formalism, point out and circumvent 
an invalidating assumption, and discuss the utility of the new, valid, relations. 

We have seen that each of the form factors F for B^ — > D^*'lv can be expanded in a 
series 

1 oo 

P(z)<P(z) 

where, in order to compare with |TjJ, we have set to — t- (N = 1). From Eq. ( |3.9| ), the 
coefficients a n obey 

al + a\ + al<\. (7.2) 

Expressing a , a±, and a 2 in terms of ^-derivatives of P(z)<f)(z)F(z) at z — 0, Eq. (|7.2|) gives 
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[P(z)cf>(z)F(z)f 



2=0 



+ 



dz 



(P(z)cf>(z)F(z)) 



2=0 



1 2 



(P(^)F(z)) 



< 1. 



2 = 



(7.3) 



For t = t_, z = corresponds to to = 1. Then Eq. ( |7.3| ) constrains the slope — p 2 and 
curvature c defined by 



F(w) = F(l) - p 2 (w - 1) + c (to - l) 2 + 0[(w- 1) 



(7.4) 



to lie within an ellipse. 

Equation ( |7.3|) is the starting point of reference |16|], with one critical simplification: 
They consider the form factor f = /q/[(Mb — Md)a/MbMd(1 and argue that it does 

not receive contributions from scalar B c mesons that would generate poles in f . This allows 
them to set P(z) = 1, leading to a nearly linear relation between p 2 and c, 



c « 0.72p 2 - 0.09. 



(7.5) 



Their argument is based on the assumption that the scalar B c mesons are broad res- 
onances because they can decay into two-particle intermediate states such as B c (0~) + rj. 
Unfortunately, the very potential models they cite refute this idea. For example, Ref. 
has the mass of the scalar 2 3 Pq state, 6.700 GeV, below the two-particle threshold for an 
r] plus the 1 % B c ground state, 6.264 + 0.548 = 6.812 GeV. The scalar 2 3 P is thus 
essentially stable with respect to hadronic transitions, since transitions involving one pion 
are suppressed by isospin, two pions by parity (or phase space, in decays to a B>*), and three 
pions by phase space. All the references |28j that calculate the relevant scalar masses and 
widths agree that there are two narrow, scalar B c resonances below the BD threshold. It is 
also worth nothing that the lowest-lying charmonium scalar state Xco is narrow. 

While sub-threshold branch cuts from states containing at least a b and a c quark may 
be legitimately ignored it is well known that poles play an essential role in the shape 
of the form factor |TT] , |55| . For this reason, the slope- convexity relations derived in are 
invalid. 
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New relations can be derived by simply including the two scalar B c states in P(t). We use 
masses from j^J, which agree with other potential model determinations to better than 1%. 
It is algebraically straightforward to input physical masses and expand our parameterization 
in powers of (w — 1), 

/oH = /o(l) 

+ (+1.72ai - 0.77/o(l)) (w-1) 

+ (-1.74ai + 0.21a 2 + 0.55/o(l)) (w - l) 2 

+ (+1.41ai - 0.27a 2 - 0.38/ (l)) (w - l) 3 

+ (-1.03a! + 0.25a 2 + 0.25/ (l)) (w - l) 4 + . . . , (7.6) 

and solve for the coefficient c of (w — l) 2 in terms of the coefficient — p 2 of (w — 1). We find 

c= 1.02p 2 + 0.21a 2 -0.23/o(l). (7.7) 

This is not a very interesting relation because the unknown coefficient a 2 , which can be as 
large as ±1 (~ ±0.6 if we include the contribution of higher states and ignore heavy quark 
symmetry violation), significantly affects the slope- convexity relationship. Had we ignored 
the Blaschke factor P(t), the coefficient of a 2 would have been 0.07. 

For N = 1, our usual truncation analysis shows that a 2 can contribute at most 4% to 
fo(w). This may seem surprising, given its obvious importance in ( [7.6|) and ( |7.7|) . The 4% 
value arises from a cancellation of the ct 2 dependence among the various (w — l) n coefficients; 
note the alternating signs of the a n coefficients. The cancellation is not accidental, but 
reflects the naturalness of an expansion in z(w;N) rather than (w — 1). This effect is 
highlighted by the observation that the (w — l) 3 term in ( [7.61) can be as large as 40% at 



«j max = 1.6, indicating that /o must be expanded to rather high order in (w — 1) if percent- 
level accuracy is desired. 

The expansion ( |7.6| ) can, alternatively, be used to test heavy quark symmetry by placing 
a restriction on the slope — p 2 . If we include the contribution from the spin- related states 
a\ is restricted to —0.61 < ai < 0.59, leading to 
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- 0.26 < p 2 < 1.84 . (7.8) 

The same relation can be derived using the form factors Sq+ or Sqq- These bounds are 
somewhat weaker than those derived from Bjorken 0] and Voloshin || inequalities^, which 
restrict 0.22 < p 2 < 1.15 once 0(a s ) corrections have been included ||38|| . 



VIII. EXPERIMENTAL FITS 

While constraints obtained from unobserved form factors like /o ma Y serve as tests of 
heavy quark symmetry, they are not well-suited for fitting to data. The reason is that 
once the truncation error on a form factor is sufficiently small (a few percent), heavy quark 
symmetry violating effects of order 20-30% become the main concern. For B — > D£u, 
using the parameterization of / + avoids any dependence on heavy quark symmetry. For 
B — > D*£u, using the form factor / minimizes the dependence on heavy quark symmetry. 
This is because on the one hand, the ratio 
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M B y/M B M D *(l + w)g 



A 



1 1 



l + —(w-2) +0 — ,— . (8.1) 



mi, ml 



2m c 

is given in terms of a single || , roughly determined constant A « 300-600 MeV and known 
f39fl perturbative functions C\jC\ = 1 + 0(a s ), while on the other hand the ratio a + /g = 
— 1/2 is determined using only spin symmetry, which is expected to hold more precisely than 
full flavor-spin symmetry [40j. We use / rather than g because it is protected from 1/m 



corrections [0] at zero recoil. 

The purpose of the following fits is not to extract the best value of \V c b\, since only the 
experimental groups themselves can correctly account for efficiencies, resolutions, smearing 



4 The physics leading to these results is quite different: The Bjorken and Voloshin inequalities 
use perturbative QCD to bound exclusive form factors directly in the semileptonic region, while 
the dispersion constraints use both perturbative QCD in the pair-production region and the phe- 
nomenological mass spectrum in the unphysical region i_ < t < t+. 
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effects, etc. Rather, since we expect the approximate results from using the QCD-derived pa- 
rameterization to survive these experimental corrections, these fits may be used to motivate 
a more thorough analysis. 



Form factors for A& — > A c £u may be extracted in the near future at CDF HTI] or LEP 



Current data for B — > D^Hu decay spectra are available from CLEO ||43||, ALEPH 



4§, OPAL @, and DELPHI @; older data exists from ARGUS @. Some progress 



has been made towards measuring individual form factors in B — > D*£u |4q| . When this 
is finally accomplished, the parameterizations for individual form factors can be applied 
without recourse to heavy quark symmetry, except for the uncertainties in the value of 
F(l). 

In the meantime, one must rely on the heavy quark symmetry prediction that the 
B — > D*iv differential rate is proportional to a function T{w) that is normalized to 
1 + 0(l/m 2 ) at zero recoil and is proportional to the Isgur-Wise function = Ef(w)f(w) 
up to l/m c corrections. We may then use the one-coefficient, QCD-derived parameterization 
of / obtained from ( |5.9| ) and to extract from data the values of | V^lJ-'(l) and a\/ 

m ~ (l + W )P f (w)M^N) t 2P ' (1 ^ (1 ' N) + 7(1) v/M.M^ J ' (8 ' 2) 
where = 1.093, z(w; N) is defined in Eq. (|4.17|) , <J)f(w;N) is given by Eq. (gjg) and 



Table 1, and Pf{w) is determined from Eqs. ( |4.25| ), ( [4.16|) , and the first four vector masses 
in Table 3. 

The procedure is precisely as detailed in ||, except that we now have only one fit co- 
efficient instead of two 0. Morever, the improvements described above reduce our one- 
parameter truncation errors to no more than 3%. We fit to the experiments whose differen- 
tial distributions are easily available. A x 2 pe r degree of freedom (d.o.f.) fit using our QCD 



dispersion relation bounds (QCD fit) to CLEO data |^3[ gives 



QCD Fit 10 3 • \V cb \ = 36.9l|? , 

= 0.000igffi , (8.3) 
Linear Fit 10 3 ■ \V cb \ ^(1) = 35.li};jj , (8.4) 
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with Xmin/ d - -f- = 0-67, while a fit to ALEPH data |4| gives 

QCD Fit 10 3 ■ \V cb \ = 31.9+H , 

yjyj = 0.096i5SS , (8-5) 
Linear Fit 10 3 ■ \V cb \ = 31.9ijj , (8.6) 

with x^ in /d.o.f. = 0.74. We have also included the values quoted by the experimental groups 
using a linear fit for comparison, and only statistical errors at one standard deviation are 
listed. 

For B — ► D£u, we can fit the parameterization of /+, 

T^y = p,»w l PM<1 ' N) + m Wtt,; N) - N)] ] • (a7) 

directly to data, without the need to invoke heavy quark symmetry. Here N = 1.108 and 
Pf + depends on the first three vector masses in Table 3. 
In this fit to ALEPH data @ yields 

QCD Fit 10 3 ■ |K6|^d(1) = 29. 2tH , 

'" 0.09518;$. 



Linear Fit 10 3 • (V^l T D (l) = 27.8±|;| , (8.9) 

withQ xLn/d-o.f. = 1.94, and now F D (1) = S /+ (l)/+(l) = 2 V ^/ + (1)/(1 + r). 

For CLEO Z? — ► D*^z/ data, which exhibits discernible curvature, our central values of 
|V^fe|jF(l) lie at the upper la boundary of the linear-fit result, while for ALEPH data, which 
is extremely flat, the central values are nearly the same. The statistical errors are larger 
because the QCD parameterization allows for curvature. 

The ALEPH B — > Div data presents an interesting area in which to test our results: 
The shape of their data will have to change as the statistical errors are reduced, if it is to 
be consistent with QCD. 



5 The larger x 2 /d.o.f. is due to their binned data point at w = 
upturn of the form factor near w max . 
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1.55, which suggests a peculiar 



IX. CONCLUSIONS 



Form factors can be reliably bounded in the pair-production region of momentum-space 
by perturbative QCD calculations. Analyticity, crossing symmetry, and dispersion relations 
may then be used to translate these bounds into constraints in the phenomenologically 
interesting semileptonic region. 

While these constraints typically imply rather weak bounds on the slopes of form factors 
njI7,n|36], quite stringent bounds can be obtained if the form factor at two or more 



points is known. The constraints actually imply an infinite number of increasingly stringent 
bounds, depending on the number of points at which the form factor is known |10|,[12| . All 
of these bounds are automatically obeyed if the form factor is parameterized as in Eq. Q3.8|) , 
even if the parameterization is truncated after a few terms. 

In this paper, we have eliminated some of the uncertainty involved in the derivation of 
these parameterizations by including two-loop perturbative corrections to the partonic side 
of the dispersion relation. We have also presented parameterizations for form factors whose 
contribution to semileptonic decay rates is suppressed by the lepton mass. 

We examined a relation between the slope and curvature of the Isgur-Wise function 
derived by Caprini and Neubert ||16|| , and pointed out the presence of sub-threshold singu- 
larities that invalidate their analysis. Once these singularities are correctly accounted for, 
the slope-convexity relations become rather weak. Bounds on the slope of the Isgur-Wise 
function made by ignoring finite mass corrections are also fairly weak. We point out that 
even with strong constraints on the slope and curvature, higher-order terms in a (w — 1) 
expansion of the Isgur-Wise function can be quite large. The parameterization in Eq. ( p.8| ) 
does not suffer from this limitation. 

For b — > c transitions, we reduced the truncation errors that describe the accuracy of 
such parameterizations. This was accomplished in part by using a parameterization that 
is automatically normalized at zero recoil to a quantity ^(1), which must be supplied by 
some other method such as heavy quark symmetry. More importantly, we included the 
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contributions of higher states in the hadronic side of the dispersion relation. These states, 
the B*D and B*D* pairs, are related to B — > D£u and B — > _D*£z/ form factors in the 
semileptonic region by heavy quark symmetry. Even in the presence of substantial heavy 
quark symmetry violation, these relations place lower bounds on the contribution from the 
higher states to the dispersion relation that lead to tighter upper bounds on the magnitudes 
of the unknown parameterization coefficients. This in turn reduces the truncation errors on 
the parameterizations of the various form factors. 

For most of the B — > D,D* form factors, the inclusion of higher states reduces truncation 
errors by roughly 40%. In one case, T\, the truncation error is reduced by as much as a factor 
of ten. This hefty improvement arises because T\ contributes very little to the dispersion 
relation, so it is far from saturating the perturbative bound until its spin-symmetry partners 
are included as well. After including all DM states, we find that each of the six form 
factors governing B — > Dlv and B — > D*£u is described to better than 3% accuracy 
using only one unknown parameter. This should be a considerable aid in experimentally 
disentangling the various form factors in differential decay distributions. 

For A b — > A c £z/ decays, there are no spin-symmetry partners to help saturate the disper- 
sion relation bound. However, the presence of more than one helicity amplitude in the same 
bound achieves the same effect. This is most dramatic in the case of Ha, whose truncation 
error is reduced by a factor of four. Of greater interest is the observation that the form factor 
F can be described at the 6% level using only one unknown coefficient. Because the \ jm c 
corrections to heavy quark symmetry relations among the baryonic form factors are given 
in terms of one additional parameter, Aa, the entire decay distribution can be described 
using only two unknown constants. This should allow a relatively clean extraction of the 
phenomenologically interesting quantity Aa- 

Finally, we used a parameterization of the form factor / and heavy quark relations to 
extract \V c b\ from B — > D*£v data. This choice of form factor minimizes the dependence on 
heavy quark symmetry. We expect the qualitative features of this extraction to persist even 
after the effects of experimental resolution, smearing, etc., are properly incorporated. Our 
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analysis suggests that the implicit error associated with the choice of parameterization is 
comparable to the statistical errors normally quoted. Similar statements apply to B — > D£u, 
except that in this case no reliance on heavy quark symmetry is necessary, since we can 
parameterize the form factor /+ directly. 

Further improvements may be possible by including additional higher states in the dis- 
persion bound, or perhaps by weighing the dispersion integral differently. One could also 
readily incorporate approximate SU(3) symmetry by using an effective ri[ = 2.5 isospin 



factor |16| in B — > D, D* decays, which would decrease truncation errors by an additional 
10%. While such improvements would be welcome for the baryonic form factors, their util- 
ity for B — > D, D* form factors is not clear. This is because, once truncation errors are 
at the few-percent level, the overwhelming source of uncertainly comes from heavy quark 
symmetry violations, which are expected to be of order 30%. Such uncertainties highlight 
the importance of extracting individual form factors, which can be parameterized using one 
coefficient without recourse to heavy quark symmetry. 



Note Added 

The essential role of the scalar B c poles in the slope-convexity relations of Caprini and 
Neubert has been pointed out independently by L. Lellouch (private communication). A 
corrected slope- convexity relation and other topics related to those in this work are in 
preparation by these authors ||49||. Additional criticism of the neglect of the scalar poles 



appears in Refs. [50 
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